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Inclusive semileptonic B decays 
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The present status of the measurement of the semileptonic B decays is reviewed. Emphasis is given on the 
factors affecting experimental errors of the Cabibbo-Kobayashi-Masawa matrix element \V c b\ and 
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1. Introduction 

In the framework of the Standard Model, the quark 
sector is characterised by a rich pattern of flavour- 
changing transitions, described by the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix. The quark tran- 
sitions B — ► XJLv and B — > X u iv provide a way for 
determining these CKM matrix elements. This re- 
port reviews the experimental measurements used to 
extract the quark mixing parameters \V c b\ and \V u b\ 
from the inclusive branching fractions of semileptonic 
b-hadron decays Br(B — > Xtv). At present, the val- 
ues of \Vqb\ obtained from inclusive semileptonics de- 
cays are more precise than the corresponding exclusive 
determinations. 

All theoretical approaches used to extract \V q b\ ex- 
ploit the fact that the mass of the b quark is large 
compared to the scale Aqcd that determines low- 
energy hadronic physics. The basis for precise cal- 
culations is the expansion in powers of Aqcd/to;, , 
where effective-field-theory methods are used to sep- 
arate non-perturbative from perturbative contribu- 
tions. 

Several studies have shown that the spectator model 
decay rate is the leading term in a well-defined ex- 
pansion controlled by the parameter Agco/mb- Non- 
perturbative corrections to this leading approximation 
arise only to order 1/m^. The key issue in this ap- 
proach is the ability to separate non-perturbative cor- 
rections, that can be expressed as a series in powers 
of l/nib, and perturbative corrections, expressed in 
owers of ct s . There are different methods Q || Q 
[j| @ to handle the energy scale fi used to sep- 
arate long-distance from short-distance physics and 
they are discussed in the theory part of these pro- 
ceedings. The coefficients of the 1/rrib power terms 
are expectation values of operators that include non- 
perturbative physics. 

In this framework, non-perturbative corrections are 
parameterized by quark masses and matrix elements 
of higher dimensional operators that at present are 
poorly known. The experimental accuracy already 
achieved, and the one expected from the large data 
sets recorded by the B-factories, make the ensuing 
theory uncertainty a major limiting factor. The ex- 
traction of the non-perturbative parameters, describ- 
ing the heavy quark masses, kinetic energy of the b 
corrections, directly from the 



quark and the I/ttk 



data has therefore become a key issue. Recently, a 
lot of progress has been made on the use of the mo- 
ments of energy and mass spectra in B — > X c lv and 
B — » A s 7 for performing these determinations. 



2. B -> XJLv 

The determination of \V c b\ from inclusive decays 
currently has an error of 2%, dominated by the 
knowledge of higher order perturbative and non- 
perturbative corrections. This method is based on 
the measurement of the total semileptonic decay rate, 
together with the leptonic energy and the hadronic in- 
variant mass spectra of inclusive semileptonic decays 
and B —> sj transitions. 

The non-calculable non-perturbative quantities are 
parametrised in terms of expectation values of 
hadronic matrix elements, which can be related to 
the shape (moments) of inclusive decay spectra. The 
shape of the lepton energy spectrum and of the 
hadronic mass spectrum provides constraints on the 
heavy quark expansion based on local Operator Prod- 
uct Expansion (OPE) that describes the properties of 
the Br(B — > X c tv) transitions. 

So far, measurements of the hadronic mass distri- 
bution and the leptonic spectrum have been made by 
Babar, Belle, Cleo and DELPHI, while CDF provide 
only the measurement of the hadronic mass spectum 

H El 

The measured hadronic mass distribution and lep- 
ton energy spectrum are affected by detector resolu- 
tion, accessible phase space, radiative corrections, etc. 
It is particularly important for reducing both theoret- 
ical and experimental uncertainties to have the largest 
accessible phase space in which to measure the decay 
spectra. Each of the respective experiments has fo- 
cused on lowering the lepton energy cut-off values. 

All experiments performed measurements of the 
hadronic invariant mass spectrum Mj^. Babar has 
measured up to the second moment of this distribu- 
tion, DELPHI up to the third. In order to compare 
with theoretical predictions, the moments are mea- 
sured with a well defined cut on the lepton momentum 
in the B rest frame. 

Belle : 8], BabarQ and CLEO 14] explored the mo- 
ments of the hadronic mass spectrum, as a function of 
the lepton momentum cuts. CLEO performs a fit for 
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the contributions of signal and backgrounds to the full 
three-dimensional differential decay rate distribution 
as a function of the reconstructed quantities q 2 , M\ , 
cos#£. Belief and Babar0 use a sample where one 
of the the B-mesons is fully reconstructed and the sig- 
nal side is tagged by a well identified lepton. In this 
case the main sources of systematic errors are the un- 
certainties related the detector and signal modelling, 
and to reconstruction. 

The hadronic mass spectrum in B — * X c lv decays 
can be split into three contributions corresponding to 
D, D* , and D**, where D** stands for any neutral 
charmed state, resonant or not, other than D and 
D*. Belle reconstructs the full hadronic mass spec- 
trum and derives M.\ with an unfolding technique. 
Babar extract the moments from the measured distri- 
butions using a calibration curve derived from Monte 
Carlo data. Employing this method Babar reaches a 
minimum momentum for the electron in the B meson 
rest frame of 0.9 GeV, while Belle reaches 0.7 GeV, 
Figure ^ 

DELPHI follow a different approach in extracting 
the moments. The D and D* component of the mass 
spectrum is well known, hence DELPHI measures 
the invariant mass distribution of the D** component 
only and fix the D and D* components. In fact, the 
D** component is not well known. For this reason 
DELPHI measures the first moment with respect to 
the spin averaged mass of D and D*. At LEP b- 
quarks were created with an energy of approximately 
30 GeV. The large boost allowed DELPHI to mea- 
sure the hadronic mass moments without a cut on the 
lepton energy [TT| . 

The shape of the lepton spectrum provides further 
constraints on the OPE. Even if these measurements 
are less sensitive to OPE parameters they are con- 
siderably more precise experimentally. Moments of 
the lepton momentum with a cut pi > 1.0 GeV/c 
have been measured by the CLEO collaboration [l5| . 
Babar ^(J extract up to the third moment of this 
distribution, using a low momentum cut of pi > 0.6 
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Figure 2: Belle:electron energy spectra for B° and 



GeV/c. Both Babar and CLEO use dilepton samples. 

Belle |3| determines the first, second and third mo- 
ments of the electron energy spectrum with a mini- 
mum lepton momentum cut ranging between 0.4 and 
1.5 GeV in the B meson rest frame. Belle measures, 
independently, the electron energy moments for the 
semileptonic decays of the B + and the B° mesons Q , 
Figure [21 In Belle, events are selected by fully re- 
constructing one of the B mesons, produced in pairs 
from T(4S), in several hadronic decay modes. Prompt 
semileptonic decays (b — > xiv) of the non-tag side 
B mesons are separated from cascade charm decays 
(B — > X c — ► Xiv), based on the correlation between 
the flavour of the tagged B and the lepton charge. 
DELPHI calculates up to the third moment for the 
lepton energy spectrum, without an explicit lepton 
momentum cut [ill ]. 

All the lepton moment measurements are consistent 
with theory and with the moment of the hadronic and 
b — ► S7 spectrum. Hadronic and lepton energy mea- 
surements are consistent within their errors. When 
compared with theory there is no sign of inconsisten- 
cies. 

A global fit to all available hadron mass, electron 
energy and photon energy moments in the IS scheme 
has been performed in [16j , giving: 



\Vcb\ = (41.4 ±0.6 ±0.1) x 10" 



(1) 




Mi(GeV-/c 4 ) 



Figure 1: Belle hadronic invariant mass spectrum 



= 4.68±0.03Ge^ 



A} s = 0.27±0.04GeF 



(2) 



(3) 



where the first error includes experimental and the- 
oretical uncertainties and the second error on \V c b\ is 
due to the B lifetime. 

A more updated set of data has been fitted in the 
kinetic scheme, resulting in [l8| : 



\Vcb\ = (42.0 ±0.2 ±0.1) x 10" 



ml™ = 4.59 ± OMGeV 



(4) 



(5) 
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Hi = 0.406 ± 0M2GeV (6) 

where the first error includes statistical and theo- 
retical uncertainties and the second error on |V^,| is 
from the estimated accuracy of the HQE for the total 
scmileptonic rate. 

In both fits, the value of the b and c quark masses 
are in good agreement with independent determina- 
tion of the same parameters. 

3. B -> XJv 

The decay rate for B — > X u lv is proportional to 
\V u t,\ 2 and m\. The theoretical description of inclusive 
B — > X u lv decays is based on the Heavy Quark Ex- 
pansion, as for B — > AT c ^i/ decays [2l],[22|, which pre- 
dicts the total decay rate with uncertainties of about 
5%. 

Experimentally, the principal challenge is to sepa- 
rate the signal B — > X u tv decays from the 50 times 
larger B — > X c iv background. This can be achieved 
by selecting regions of phase space in which this back- 
ground is highly suppressed. In these regions the 
spectra are affected by the distribution of the fa- 
quark momentum inside the B meson, which can be 
described by a structure or "shape function" (SF) 
[23l |24| , in addition to weak annihilation and other 
non-perturbative effects. 

Extrapolation from the limited momentum range 
near the endpoint to the full spectrum is a difficult 
task. The shape function is a universal property of B 
mesons at lea ding order |23l l24j . Sub- leading shape 
functions H^|-j33| arise at each order in 1/rrib and 
differ between semileptonic and radiative B decays. 
Several functional forms for the SF, which generally 
depend on two parameters related to the mass and 
kinetic energy of the b-quark, A or nib > and Ai or 
7r 2 , have been proposed. The values and precise def- 
initions of these parameters depend on the specific 
ansatz for the SF, the mass renormalization scheme, 
and the renormalization scale chosen. 

In inclusive measurements, the most common kine- 
matic variables discussed in the literature, each having 
their own advantages, are the lepton energy (E e ), the 
hadronic invariant mass (Mx), the leptonic invariant 
mass squared (q 2 ) and the light-cone momentum com- 
ponent P+ = Ex — \ Px\- In all cases, the experiments 
need to model B — » X u lv decays in order to calculate 
acceptances and efficiencies. 

3.1. Endpoint measurements 

The first \V u b\ measurement performed with this 
technique was by CLEO. In the rest frame of the B 
meson, the kinematic endpoint of the electron spec- 
trum is about 2.3 GeV/c for the dominant B — > X c lv 
decays and about 2.6 GeV/c for B — > X u lv decays. 




1.1 1.5 1.9 2.3 2.7 3.1 3.5 
Electron Momentum (GeV/c) 

Figure 3: BaBar: Electron momentum spectra in the 
T(4iS) rest frame: (a) on-resonance data (open circles) and 
scaled off-resonance data (solid circles); (b) on-resonance 
data after non-BB background subtraction (triangles) 
compared to simulated BB background (histogram); (c) 
on-resonance data after subtraction of all backgrounds 
(data points), compared to the simulated B — > X u lv sig- 
nal spectrum (histogram). 



The primary challenge in reducing the lepton momen- 
tum cut in the endpoint method is controlling the 
B — > X c lv background at the required precision. The 
spectrum above 2.3 GeV/c is dominated by electrons 
from B — > X u iv transitions, and this allows for a 
relatively precise measurement, largely free from BB 
background, in a 300 MeV/c interval that covers ap- 
proximately 10% of the total electron spectrum for 
charmless semileptonic B decays. Figure shows the 
Babar lepton energy spectrum distribution for the sig- 
nal and background. Belle extracts the B — > X u lv sig- 
nal in the momentum region 1.9 — 2.6 GeV/c, Babar 
covers 2 - 2.6 GeV/c and CLEO 2.3 - 2.6 GeV/c 

mmip. 

3.2. M x , q 2 and P+ 

As described earlier, this method relies on a sam- 
ple of events where one of the B mesons from the 
T(45) decay is fully reconstructed, while the semilep- 
tonic decay of the signal side B meson, is identified by 
the presence of a high momentum electron or muon. 
With this sample we are able to construct all kine- 
matic variables: the invariant meson mass Mx, the 
lepton neutrino mass squared q 2 , and the hadronic 
light-cone momentum P + , which are the best avail- 
able discriminators of signal and background in inclu- 
sive \V u b\ analyses. 

BABAR selected semileptonic B — > X u lv decays 
using an approach based on simultaneous require- 
ments for the electron energy, E e , and the invari- 
ant mass squared of the ev pair, q 2 [35| . The neu- 
trino 4-momentum is reconstructed from the visible 
4-momentum and knowledge of the e + e~ initial state. 
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Figure 4: The electron momentum spectrum in the T(45) 
restframe: (a) P+ distribution for the selected events, (d) 
P+ distribution after subtracting the B —* X a tv contribu- 
tion. 



The dominant charm background is suppressed by se- 
lecting a region of the q 2 -E e phase space where cor- 
rectly reconstructed B — > X c lv events are kinemati- 
cally excluded. 

Theoretical studies indicate that it is possible to 
reduce the theoretical error on the extrapolation by 
applying simultaneous cuts on Mj^and q 2 in inclusive 
B — ► X u lv decays [3(|. In fact, while the Mx dis- 
tribution has a large usable fraction of events, of the 
order of 70%, but depends on the shape function de- 
scribing the Fermi motion of the b quark inside the 
B meson, the q 2 distribution is less sensitive to non- 
perturbative effects and less dependent on the calcu- 
lation. Unfortunately, only a small fraction of events 
(about 20%) are usable with a pure q 2 selection. The 
study in [36| shows that combined cuts on Mx and 
q 2 mitigate the drawbacks of the two methods while 
retaining good statistical and systematic sensitivities. 
Babar performed a measurement of the \V u b\ CKM 
matrix element, on the fully reconstructed B sample, 
by using the combined information of the Mx-q 2 dis- 
tribution to discriminate signal and background and 
to minimize the theoretical uncertainties 37]. 

Belle has measured partial rates with cuts on Mx, 
q 2 , E e , but the highlight is the first measurement 
of the light-cone momentum P + |38( ■ Belle measures 
the partial branching fraction in the kinematic region: 
P+ <0.66 GeV/c, Figure 0] The direct measurement 
of P+ is quite important as the shape function be- 
comes influential when P + is not large compared to 



3.3. Reduced Model Dependence 

Babar used two new techniques to extract | V u b | from 
inclusive semileptonic B decays where the uncertain- 
tics due to mb and the modeling of the Fermi motion 
of the b quark inside the B meson are significantly 
reduced [39 ]. 

Leibovich, Low, and Rothstein (LLR) presented 



a prescription to extract |Ki&| with reduced model 
dependence from either the lepton energy or the 
hadronic mass Mx [4(j ■ A technique to utilise weight 
functions had been proposed previously by Neubert 
|4l| The calculations of LLR are accurate up to cor- 
rections of order a 2 and (Ams/(£m;,)) 2 , where £ 
is the experimental maximum hadronic mass up to 
which the B — > X u lv decay rate is determined and 
A ~ AQCD. This method combines the hadronic 
mass spectrum, integrated below £ , with the high- 
energy end of the measured differential B — > sj pho- 
ton energy spectrum via the calculations of LLR. An 
alternative method ^2 to reduce the model depen- 
dence is to measure the B — > X u lv rate over the en- 
tire Mx spectrum. Since no extrapolation is necessary 
to obtain the full rate, systematic uncertainties from 
m,b and Fermi motion are considerably reduced. Per- 
turbative corrections are known to order a 2 . Babar 
extracts the B — * X u h> rate from the hadronic mass 
spectrum up to £ = 2.5 GeV/c2 which corresponds to 
about 96% of the simulated hadronic mass spectrum, 
Figure |SJ 
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Figure 5: \V u b\ as a function of £ with LLR method (left) 
and for the determination with the full rate measure- 
ment (right). The errors are statistical only. The shaded 
area illustrates the theoretical uncertainty; the inner light 
shaded (yellow) area indicates the perturbative part of the 
uncertainty. The arrow indicates £ = 1.67GeV/c 2 . 



Using the weighting technique with the photon en- 
ergy spectrum in B — > in Ref. 43] the hadronic 
mass spectrum up to a value of £= 1.67GeV/c2, Babar 
finds \V ub \/\V ts \ = 0.107 ± 0.009 s t at ± 0.006 syst ± 

0. 007th eo - Assuming the CKM matrix is unitary gives 
\V U \ = \V cb \ x (1 ± 0(1%)) and taking \V cb \ from 
one gets: \V ub \ = (4.43 ± 0.38 ± 0.25 ± 0.29) x 10" 3 , 
where the first error is the statistical uncertainty, the 
second is systematic, the third theoretical. Babar 
also determines | V^t, | from the full Mx spectrum, 

1. e., up to a value of £=2.5GeV/c 2 , and find \V llb \ = 
(3.84 ± 0.70 sfat ± 0.30 sys ± 0.10 tfeeo ) x lO" 3 . 
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Figure 6: Kinetic mass fit from which the heavy quark pa- 
rameters needed in calculating B —> X u lv partial rates are 
derived. The dotted curve shows the fit using information 
form B — > X c lv only, the red line includes information 
from B — > s'y. 



3.4. HQE parameters and shape function 
input 

The global fits to B — > XJLv and B — > S7 mo- 
ments in ref discussed earlier provide input values 
for the heavy quark parameters needed in calculating 
B — v X u lv partial rates and to constrain the first and 
second moments of the shape function. Additional 
information on the leading shape function and HQE 
parameters is obtained from the photon energy spec- 
trum in B — > S7 decays. The results of the global fit 
are shown in Figure El 

3.5. extraction 

The CKM parameter \V u b\ is obtained directly from 
the partial branching fraction using 

\V ub \ 2 = AB(B -> X u £p)/(Rt b ), 

where tb is the average B lifetime. R is the theoretical 
prediction of the partial rate and is calculated for a 
given signal region by using an inclusive B — ► X u iv 
decay generator. The extracted values of \V u b\ are 
given in Figure 

The world average is determined by HFAG. We have 
chosen the value extracted using the BLNP j^j the- 
oretical framework as the world average of the CKM 
parameter \V u b\. The value of m& used in the world 
average is from the global fit in the Kinetic scheme to 
measurements of B — > X c £u and B — > decays. 

HFAG also extracts \V u b\ using the Dressed Gluon 
Exponentiation |45j| framework, a recent new addition 
to the phenomenology landscape of inclusive B-meson 
decays. In this framework the on-shell b-quark cal- 
culation, converted into hadronic variables, can be 



Figure 7: The world average for \V u t\ 



directly used as an approximation to the meson de- 
cay spectrum, without need of a leading-power non- 
perturbative function i.e no shape function. The on- 
shell mass of the b-quark within the B-meson is re- 
quired as an input (m&). Theoretical uncertainties 
are assessed by varying nib, the strong coupling con- 
stant (a s ), the number of light fermion flavours, and 
the method and scale of the matching scheme intrin- 
sic to the approach. As the DGE error is not yet fully 
calculated (e.g. the weak anihilation contribution is 
missing), we prefer not to quote this method as the 
principal one. It must be noted that the value of \V u b\ 
obtained within the DGE framework agrees very well 
with the BLNP determination. 

At present, as indicated by the average given above, 
the uncertainty on \V u b\ is at the 7% level. The un- 
certainty on rrib used here is 40 MeV, contributing an 
uncertainty of 4.5% on \V u b\- 



4. Conclusions 

The global fits to the semileptonic and radiative 
decay spectral moments are used to extract the CKM 
matrix element \V c b\ with 2% precision dominated by 
HQE measurements. 

Various methods to extract \V u b\ are used, giving 
a 7% precision in the BLNP scheme, dominated by 
HQE parameters. 
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